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ABSTRACT 

We present a mark correlation analysis of the galaxies in the Sloan Digital Sky Survey using weights 
provided by MOPED. The large size of the sample permits statistically significant statements about 
how galaxies with different metallicities and star formation histories are spatially correlated. Massive 
objects formed a larger fraction of their stars at higher redshifts and over shorter timescales than did 
less massive objects (sometimes called down-sizing). We find that those galaxies which dominated the 
cosmic star formation at z « 3 are predominantly in clusters today, whereas galaxies which dominate 
the star formation at z w inhabit substantially lower mass objects in less dense regions today. Hence, 
our results indicate that star formation and chemical enrichment occured first in the denser regions 
of the Universe, and moved to less dense regions at later times. 

Subject headings: stellar populations - large scale structure of universe 



1. INTRODUCTION 

There has been substantial recent progress in the de- 
velopment of methods which determine the star forma- 
tion and chemical composition histories of galaxies from 
the integrated spectra of their stellar populations. The 
traditional approach has been to determine the instan- 
taneous star formation rate or metallicity from certain 
features in the spectrum of a galaxy. However, sev- 
eral recent algorithms ( Heavens et all 200(]HSodre et al.1 



l2TO5t iMathis et al.ll2006HOcvirk et al.ll2006|) . have been 
developed which use the entire spectrum to infer the 
entire star formation history and the evolution of the 
chemical composition of the object, rather than simply 
the instantaneous values of these quan tities. One such 
method, MOPED l|Heavens et alJl200(]|) . has been used 
to determine the star formation histories and metallici- 
ties of galaxies drawn from the Sloan Digit al Sky Survey 
jPanter et al." 2003; Heavens et al. 2004; "Pante r et all 
|2D04; Jimenez ct al. 2005; Panter et al. 2006a b). 

There has also been significant progress in quantify- 
ing how galaxies are distributed on large scal es, and us- 
ing this to constrain cos mological parameters IjCole et al.l 
120051: ISeliak et al]l2005D . In such analyses of galaxy clus- 
tering, it is common to treat galaxies as points, ignoring 
the fact that galaxies have different luminosities, colors, 
masses, star formation histories, metallicities, etc. How- 
ever, as a result of improvements in detector technol- 
ogy, and in the algorithms such as MOPED with which 
the new data is analyzed, many such galaxy attributes 
are now sufficiently reliably measured that one can use 
them as weights when studying the clustering of galax- 
ies. Thus, one can now study the clustering of lumi- 
nosity, color, star formation rate, etc. Mark statistics 
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IjStovan fc Stovanlll994l: iBeisbart fc Kerscheill2000D pro- 
vide a useful framework for describing attribute- weighted 
clustering. Moreover, they provide sensitive probes of 
how the pr operties of g alaxies correlate with their envi- 
ronments (Sheth "200^. In this respect, mark statistics 
provide a useful link between the large-scale structures 
which galaxies trace, and the properties of those galaxies. 
They have recently been used to measur e the clusterin g 
of luminosity and color in the SDSS (Skib ba et al.l20 0fil. 
In this Letter, we use mark correlations to demonstrate 
the wealth of information which the MOPED-determined 
histories provide. We show that a mark correlation anal- 
ysis, using MOPED attributes as marks, allows one to 
study how the star formation histories of objects are cor- 
related with their present-day environments. Our main 
finding is that the majority of close (< 2 Mpc) galaxy 
pairs today are made up of objects which formed the 
largest fraction of their stars at z « 3, whereas those 
which formed stars at an above average rate within the 
last Gyr or so tend to be in less clustered environments 
today. In hierarchical models, overdense regions today 
were overdense in the past. Thus, our results indicate 
that star formation at high-redshift occured in dense re- 
gions, whereas it only occurs in substantially less dense 
regions today. A similar study of the metallicity shows 
that the close pairs which formed a larger than aver- 
age fraction of their stars at z « 3 also have above av- 
erage metallicities, but that there is no correlation be- 
tween metallicity and environment for stars formed more 
recently. Where necessary, we assume a flat ACDM 
model with fjn = 0.27 and Hq = 71 km s'^Mpc-^ 
llSpergelet al.ll200l . 

2. MEASUREMENTS IN THE SDSS 

As described in iPanter eraTI l|2006albD . the MOPED 
algorithm has been used to extract star formation 
and metallicity histories of a magnitude limited sample 
(15.0 < rur < 17.77) of about 300,000 galaxies drawn 
from the Third Data Release of the Sloan Digital Sky 
Survey (SDSS DT^,3: lAbazaiian et a,].l (12003) 'I. We mea- 
sured mark correlations, using MOPED attributes as the 
marks, in two smaller volume-limited catalogs extracted 
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Fig. 1. — Mark correlations in the two volume limited catalogs 
described in the text; symbols with error bars show results for 
the fainter sample, and lines without error bars are for the more 
luminous sample. Close pairs tend to be more luminous (top panel) 
and to host older stellar populations (bottom) than average. 



from the SDSS DR3. To facilitate future halo- model 
based interpretations of our measurements, these cata- 
logs were chosen to app roximately correspon d to those 
studied bv lZehavi et alJ l|2005() and lSkibba ~ al. (2006): 
a brighter catalog which spans Mj. < —21.5 for which 
0.02 < z < 0.135, and one which includes fainter ob- 
jects, Mr < —20, and so spans a smaller redshift range: 
0.02 < z < 0.071 (the associated apparent magnitude 
limit is conservative: rur < 17.5). 

Figure Q shows the mark correlations in the two cat- 
alogs when r— band luminosity (top) and the MOPED- 
inferred mass-weighted age (bottom) arc the marks (sym- 
bols with error bars show results for the fainter sample, 
and lines without error bars are for the brighter sam- 
ple). So that the weights are dimensionless, the weight 
of each galaxy is normalized by the mean value for the 
population. The notation WW/DD indicates that the 
mark statistic is the ratio of weighted pair counts to un- 
weighted pair counts. (In this notation, the traditional 
unweighted correlation function would be DD/RR, where 
RR is the number of unweighted pair counts in a ran- 
dom distribution.) Symbols with error bars show results 
for the fainter catalog, and lines without error bars are 
for the more luminous sample — we use a similar conven- 
tion in all the figures which follow. Error ba rs were esti- 
mated using t he analytic expression s given in lSheth et al.l 
l|200(it) . which 'Skib ba et aLl 12006") have shown are sim- 
ilar to those from a jackknife analysis; they are similar 
for the two samples. The results of Skibba et al. also 
show that, while using redshift- rather than real-space 
distances tends to make WW/DD closer to unity on small 
scales, this is not a severe effect. 

The top panel shows that close pairs of galaxies tend to 
be more luminous than average, consistent with previous 
mark correlation analyses of SDSS galaxies l)Skibba et al.l 
12006(1 . Halo- model based analyses of the clustering 
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Fig. 2. — Joint distribution of star-formation fraction and metal- 
licity in the volume limited catalog with Mr < —20; Different 
panels show results for different lookback times, t in Gyrs, and as- 
sociated redshifts, z. In all panels, the marks have been normalized 
by the mean mark in the bin. 



of SDSS galaxies indicate that, on scales smaller than 
1 Mpc, the pair c ounts are dorn inated by galaxies in 
massive haloes (Z ehavi et al.ll2005|) . Thus, the top panel 
shows that galaxies in massive haloes are over-luminous. 

The bottom panel shows that, in addition to being 
more luminous than average, close pairs tend to have 
older than average stellar populations. This suggests 
that the most massive halos host the oldest stellar popu- 
lations. The scale dependence in the bottom panel is also 
qualitatively simila r to tha.t seen in semi-analytic g alaxy 
formation models IjShethl l2005t iSheth et alJl2006tl . In 
the models, this happens because close pairs are domi- 
nated by galaxies in clusters, and cluster galaxies host 
the oldest stars. This is consistent with the halo-model 
interpretation mentioned above. 

Figure n uses marks which are the result of integrating 
over the entire star formation history of each object. One 
of the great virtues of the MOPED analysis is that it 
returns not just the mass-weighted age at the present 
time, but an estimate of the entire star formation history 
of an object. Thus, for each object, we have constructed 
estimates of the fraction of the current stellar mass which 
formed in each of eight bins in lookback time. 

Figure |2] shows the joint distribution of star forma- 
tion fraction and metallicity in the fainter catalog, for 
each of the lookback time bins. In each panel, the 
marks have been normalized by the mean value in the 
bin. For instance, at the two largest lookback times, 
((SFF),(Z/Z0)) = (0.25,0.65) and (0.43,0.77). These 
numbers are (0.25,0.70) and (0.55,0.89) in the more lu- 
minous catalog. The differences indicate that massive ob- 
jects formed a larger fraction of their stars at higher red- 
shifts, on a shorter t imescale, than lo wer mass objects — a 
point made bv .Heavens et alJ 1(20041) . 

In all panels there is a population of objects which 
have small star formation fractions but a large range of 
metallicities, and a population which has small metal- 
licities but a tail which extends to large star formation 
fractions. However, at large lookback times, there is an 
additional population which has above average metal- 
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Fig. 3. — Mark correlation functions, with star-formation fraction as the mark, in the two volume limited catalogs described in the text. 
Symbols with error bars show results for the fainter sample; lines without error bars represent the more luminous sample. Panels show 
results for different lookback times, t in Gyrs, with associated formation redshifts, z. In both catalogs: close pairs today had higher than 
average star formation fractions at z = 2.5, average star formation fractions at 2 = 1, and lower than average star formation fractions more 
recently. The anti-correlation between star formation fraction and environment persists up to lookback times of 5 Gyrs. 



licity and star formation fraction. For the two largest 
lookback times, this population comprises 24% and 39% 
of the galaxies in the fainter catalog (26% and 51% in 
the brighter catalog). Where are these objects now? 

Figure |31 shows a mark correlation analysis of the star 
formation fraction at these eight lookback times. The 
bottom right panel shows that close pairs today had 
larger than average star formation fractions 11 Gyrs ago. 
Comparison with the other panels indicates that these 
pairs had average star formation fractions at redshifts of 
order unity, and smaller than average star formation frac- 
tions more recently. Thus, our analysis provides graphic 
evidence that the objects which underwent vigorous star 
formation at the highest redshifts are currently in clus- 
ters, where the current star formation rate is smaller than 
average. 

Figure 0] shows a similar analysis when metallicity, 
Z/Zq, is the mark. The bottom right panel shows that 
the close pairs which had above average star formation 
fractions at z = 2.5 also tend to have above average 
metallicities. There are no clear correlations with en- 
vironment in the other panels, except, possibly, at the 
smallest lookback times (top left), where there is an en- 
hancement on scales where the pair counts are dominated 
by galaxy groups. Thus, our measurements indicate that 
the stellar populations of the most massive halos are old 
and metal rich. The population of objects which had 



above average star formation fractions and metallicities 
at large lookback times (c.f. Figure EJ are today in clus- 
ters. 

3. CONCLUSIONS 

A mark correlation analysis of SDSS galaxies using 
MOPED-derived ages, metallicities and star formation 
histories shows that close pairs tend to host stellar pop- 
ulations which are older than average (Fig. ^. Close 
pairs also tend to have formed a larger fraction of their 
stars at z w 3 than average, but the star forming frac- 
tion at z < 1 of such pairs is below average (Fig.|2Jl. The 
objects which were forming stars at z « 3 have above 
average metallicities (Fig. 0J . These trends do not de- 
pend significantly on the mean luminosity of the sample, 
so they are approximately independent of stellar mass. 

Close pairs are dominated by galaxies in massive halos. 
Hence, our results indicate that galaxies which formed 
the bulk of their stars at high redshift are today in clus- 
ters, in which there is little ongoing star formation. Since 
clusters formed from overdense regions in the early Uni- 
verse, our results imply that cosmic star formation has 
moved from dense to ever less dense regions. This is qual- 
itatively consistent with the findings of Poggianti et a3 
(J2006) . However, whereas Poggianti et al. measure in- 
stantaneous star formation rates in cluster galaxies iden- 
tified over a wide range of redshifts (the SDSS at 2; ^^ 0, 
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Fig. 4. — Mark correlation functions with mctallicity as the mark in the two volume limited catalogs described in the text. Symbols 
with error bars show results for the fainter sample, and lines without error bars are for the more luminous sample. Panels show results 
for different lookback times. Close pairs which formed their stars 11 Gyrs ago tend to have metallicities which are above average for that 
epoch (bottom right). At smaller lookback times there is little correlation between metallicity and present day environment except very 
recently (top left) where most of the enrichment is taking place in galaxy groups. 



and the ESO Distant Cluster Survey for 0.4 < z < 0.8), 
our method uses the spectra of the local galaxy popula- 
tion to infer the entire cosmic star formation history. In 
particular, it does not require classification of the galax- 
ies into 'cluster' and 'field' populations, nor docs it re- 
quire acquisition of a galaxy sample which spans a wide 
redshift range. It is remarkable that these two very dif- 
ferent methods agree. 

For similar reasons, the top left panel of Fig. 13 may be 
compared with recent studies of the dependence of cur- 
rent star formation on environment ( Gomez ct al. 2003; 
IRalogh et all l2fiol IMa.teus et all l200(1) . We all find 
smaller star formation rates in dense regions today. Note, 
however, that our analysis is not restricted to the current 
epoch — it covers 11 Gyrs in lookback time. Halo model 
interpretations of our measurements will help determine 



if the environment plays a crucial role in regulating star 
formation. 
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